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It has been recently discovered that stabilization of two-dimensional (2D) solitons against the 
critical collapse in media with the cubic nonlinearity by means of nonlinear lattices (NLs) is a chal- 
lenging problem. We address the ID version of the problem, i.e., the nonlinear-Schrodinger equation 
(NLSE) with the quintic or cubic-quintic (CQ) terms, the coefficient in front of which is periodically 
modulated in space. The models may be realized in optics and Bose-Einstein condensates (BECs). 
Stability diagrams for the solitons are produced by means of numerical methods and analytical 
^ , approximations. It is found that the sinusoidal NL stabilzes solitons supported by the quintic-only 

■ nonlinearity in a narrow stripe in the respective parameter plane, on the contrary to the case of 

the cubic nonlinearity in 2D, where the stabilization of solitons by smooth spatial modulations is 
not possible at all. The stability region is much broader in the ID CQ model, where higher-order 
solitons may be stable too. 
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" PACS numbers: 42.65.Tg, 05.45.Yv, 03.75.Kk, 03.75.Lm 
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^| I. INTRODUCTION 

The use ofperiodic potentials, induced by photonic crystals or optical lattices for controlling light beams in optical 
waveguides [l[ or matter waves in Bose-Einstein condensates (BECs) Q was the subject of many works, see reviews 
It has been demonstrated that lattice potentials balancing the self- focusing or defocusing cubic nonlinearity 
uphold bright solitons — ordinary ones or gap modes, respectively. The general topic of solitons in period potentials 
has been reviewed in detail in Refs. In the two-dimensional (2D) geometry, where the cubic (x*^) self- focusing 

leads to the critical collapse @, fundamental solitons in free space (Townes solitons Q) are unstable, representing a 
separatrix between collapsing and decaying modes. Vortex solitons in the 2D free space exist too [101 ]. being subject 
fNJ ■ to a still stronger instability against azimuthal perturbations [1, HH). However, lattice potentials can stabilize 2D 
' fundamental solitons and solitary vortices alike [12]. Moreover, the same potentials may also stabilize 2D solitons 
against the supercritical collapse, driven by the self- focusing quintic (x ) term |13j . 
— i ■ In ID, the critical collapse is induced by the quintic nonlinearity, which gives rise to the known ID counterparts 
of the Townes solitons, see Eq. @ below. The analogy to the 2D setting suggests a possibility of the stabilization 
of the ID variety of Townes solitons by means of the single-well or lattice potentials, which was demonstrated in 
Refs. [1-4] and pl|, respectively, with the help of numerical simulations and a variational approximation (VA). As 
. concerns physical realizations, the quintic nonlinearity (acting, generally, in a combination with the cubic term) was 
theoretically predicted fill an d experimentally observed in diverse optical media, including liquids [13, HU ; glasses 
fl9j and ferroelectric films [20j |. Generally, this nonlinearity comes along with the cubic losses. Nevertheless, the losses 
may be usually neglected over relatively short propagation distances for which experiments are run. In particular, 
the cubic-quintic (CQ, alias x^ '■ X^ 5 ') optical nonlinearity may be efficiently adjusted in colloids formed of metallic 
nanoparticles, by selecting the size of the particles and the colloidal filling factor pjj]. Using these parameters, one 
may emulate a medium with the quintic-only nonlinearity, i.e., x^ — 0, which is the case especially interesting for 
the theoretical consideration, see below. Another realization of an effectively quintic medium may be possible in a 
composite material built as an alternation of narrow self- focusing and defocusing layers [2l|, in which the cubic x'' 3 '' 
coefficient is canceled on the average, while its x^ coefficient remains uncompensated. 

In terms of the BEC, the quintic term in the Gross-Pitaevskii equation (CPE) may represent three-body interactions 
in a sufficiently dense condensate, provided that the lossy part of these interactions may be neglected (the latter 
depends on particular properties of the condensate, the use of the Feshbach resonance for controlling the scattering 
length of inter-atomic collisions, etc.) [22j . On the other hand, a universal quintic term appears in the effective 
GPE, in the absence of any three-body interactions, as a result of the reduction of the dimension from 3 to 1 under 
the action of the strong confinement in the transverse plane [23|, [24j. In the latter case, no dissipative nonlinearity 
appears, and the sign of the effective quintic term is self-attractive, irrespective of the sign of the scattering length 
which determines the self-attractive or repulsive sign of the cubic nonlinearity. The relative strength of the respective 
quintic term in comparison to the basic cubic one can be estimated as ~ 5|a|niD (see, e.g., Ref. |24j), where a is the 
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scattering length (a < corresponds to the attraction between atoms), and nio the ID density of the condensate. For 
the typical situation with matter- wave solitons containing ~ 5000 atoms in the area of size ~ 2 /im, with a ~ —0.1 
nm [25| . the estimate yields the value ~ 1, hence the quintic terms may be as important as the cubic one. 

More recently, the studies of the soliton dynamics in periodic potentials were extended for nonlinear lattices (NLs), 
which are represented by pseudopotentials [26[ induced by spatially periodic modulations of the local strength of the 
nonlinearity, see original works [27j|-[33| and recent review In optics and BEC, NLs may be engineered by means 
of various techniques, such as filling voids of photonic crystals with properly selected materials, or using the Feshbach 
resonance controlled by spatially patterned fields |7[ (in particular, magnetic lattices governing the local strength of 
the Feshbach resonance have been built using sets of ferromagnetic films [2!|). In quasi-lD BECs, an NL can also 
be induced by periodically modulating the transverse-confinement strength in the axial direction [30j. Similarly, a 
periodic modulation of the local nonlinearity strength in a planar optical waveguide may be created by varying the 
thickness of the waveguiding layer. In particular, the periodically modulated quintic-only and CQ nonlinearities, 
considered below, may be realized in the above-mentioned colloid layers with the periodic variation of the transverse 
thickness, or by applying a properly patterned magnetic field to the colloidal suspension formed by ferromagnetic 
particles. 

The previous analysis has demonstrated that NLs readily support stable ID solitons, but it is quite difficult, although 
possible, to stabilize 2D solitons against the critical collapse by means of nonlinear (pseudo)potentials. Namely, 2D 
sinusoidal lattices completely fail to provide for the stabilization, while shapes with sharp edges, such as circles or 
stripes, are able to do it [3l|-[33|. Although a periodic lattice of nonlinear circles with sharp edges can support stable 
2D solitons, they are virtually the same as those supported by a lone circle, i.e., the periodicity plays no essential role 
in the stabilization (32|. 

The results outlined above suggest to investigate the ID counterpart of the problem, i.e., stabilization of ID solitons 
by means of NLs against the respective critical collapse, accounted for by the x^ 5 ' (quintic-only) or CQ nonlinearity. 
Below, we perform the analysis of the model with the quintic-only and CQ self-focusing terms, whose strength is 
subject to the periodic spatial modulation, by means of numerical methods and analytical approximations. We find 
that a smooth sinusoidal NL can stabilize ID solitons, on the contrary to the negative result in 2D, but only in a 
narrow parametric region. This finding is explained with the help of an analytical approach. The sinusoidal CQ NL 
gives rise to a much broader stability area than its quintic counterpart, which is demonstrated by means of the VA, 
in addition to numerical results. A chain of nonlinear domains with sharp edges may support stable solitons in ID 
as well, but it can be readily checked that, like their 2D counterparts, such solitons are actually supported by lone 
nonlinear islands, the effect of the periodicity being negligible. 



II. THE MODEL AND ANALYTICAL APPROXIMATIONS 



In terms of optical realizations, the scaled form of the underlying nonlinear Schrodinger equation (NLSE) for the 
field amplitude, u{x,z), is 

iu z = -~u xx - [e + £i cos(2x)] |u| 2 u - [50 + 51 cos(2x)}\u\ i u, (1) 

where £0 and 50 are strengths of the uniform CQ terms, e\ and g± represent the NL, and the scale of coordinate 
x is fixed by setting the NL period to be n. Below, the remaining scaling invariance is used to fix g x = 1. The 
corresponding GPE for the mean-field wave function of BEC is obtained from Eq. ([T]) by replacing propagation 
distance z with time t. A model featuring the periodic modulation of the cubic nonlinearity along the evolution 
variable (z), rather than along x, was studied too, supporting stable motion of solitons in the discrete version of the 
ID model 

Stationary solutions to Eq. (JT]) with propagation constant k are sought for as u(x. z) — U(x)exp{ikz), where real 
function U(x) satisfies equation 

kU = hj xx + [e + £1 cos(2a:)] U 3 + [g + g x cos(2x)]t/ 5 , (2) 
which can be derived from the Lagrangian, 



L = — 



+ OG 



kU 2 + t:UI- -[eq + sx cos{2x)} U 4 



-3 [ffo cos(2a;)] U 6 } dx. 



(3) 
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In the case of £0,1 = .91 = 0, the above-mentioned unstable Townes-type solution of the quintic NLSE with go > is 
well known: 



U 2 (x) = v/3fc/ffosech (y/8kx) , 



whose norm (the total power or number of atoms, in terms optics or BEC, respectively) does not depend on k, 

V= / U 2 (x)dx = ^J^-, 
2 V 2g 



(4) 



(5) 



which is the characteristic feature of Townes solitons [8[ . The addition of an arbitrarily weak cubic self- focusing term, 
with constant coefficient e > 0, lifts the degeneracy and turns solution (Q} into exact stable solutions of the CQ 
equation [H, |36j]. The addition of the self-defocusing cubic term with £0 < lifts the degeneracy too, but makes the 
entire soliton family unstable [24|, . 

To apply the VA [13], we adopt the Gaussian ansatz for the soliton with width W and norm N: 



U 2 (x) 



N 



exp 



x 



The substitution of the ansatz into Lagrangian ^ yields the effective Lagrangian, 
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which gives rise to the corresponding variational equations, 8L c s/dW — dL c g/dN = 0: 
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As shown below, the VA produces reasonable results in the model with the CQ nonlincarity. The VA may be extended 
to investigate the stability of perturbed (nonstationary) solitons by including the chirp into the ansatz, but this is 
not really necessary, as the stability of the solitons which are reasonably well approximated by the ansatz is correctly 
predicted by the VK criterion, as shown below. 

The VA fails to produce accurate results for the quintic-only model (eo,i = 0), unlike the CQ one (see below). 
Nevertheless, in the quintic model another analytical approximation can be applied to broad solitons with k < 1 and 
g 3> g\. In this case, an approximate correction to soliton generated by term ~ gi in Eq. ([5]), is found to be 



(31/2) (3fc/ 5o ) 5/4 sech(V8fcx) cos (2a;), 



and the corresponding correction to norm ([5]), which lifts its degeneracy, is 



N 1 » (3V3tt/4) .gi. 9o 3/2 exp {-%/ (2\/2fc)) 



(8) 



(9) 



Thus, correction ((HJ replaces the single value ([5]) by the exponentially narrow stripe of values ([5]), within which the 
norm satisfies the Vakhitov-Kolokolov (VK) necessary stability criterion, dN/dk > [1, [38| , suggesting that there 
may exist a narrow stability stripe for the ID solitons in the quintic-NL model. This expectation is confirmed by 
numerical results presented below. 

For the sake of the comparison of the models with the critical (quintic and cubic, respectively) nonlinearities in ID 
and 2D, it is relevant to mention that an approximation similar to that based on Eq. ([8]) can be developed too for 
broad 2D solitons in the model with the full 2D or quasi-lD sinusoidal modulation of the nonlinearity, making use 
of formula J Q 27r cos (2r cos#) d6 = 2ttJq (2r) for the angular integration in the calculation of the ensuing correction for 

the norm, and of the asymptotic form for the 2D Townes soliton, Uo{r) ~ r~ 1 / 2 exp \/2krj . However, the eventual 

result for the 2D setting turns out to be opposite: dN^ 2D ^ / dk < 0, which explains the failure to of the sinusoidal NL 
to stabilize the 2D solitons. 



4 




4 




(b) 


=> 2 









--^ ■ 





-4 -2 




X 



FIG. 1: Shapes of stable (a) and unstable (b) solitons in the quintic model with e = 0, go = 0.5 are shown for k — 69.8, N = 1.65, 
and k = 71.5, N = 1.78, respectively. The dashed lines represent the shapes predicted by the variational approximation. 
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FIG. 2: (Color online) The evolution of the stable and unstable solitons from Fig. [JJis displayed in the left and central panels, 
respectively. In addition, the right panel shows an example of the unstable evolution towards the collapse in the quintic model 
(unlike the decay in the central panel), for the soliton with k = 88.2 and TV = 1.85. 



III. LOCALIZED MODES IN THE QUINTIC MODEL 



Stationary solutions to Eq. (pj, i.e., solutions of Eq. ([2]), were found by means of the imaginary-time-propagation 
method [3^ |. and their stability was tested by subsequent simulations of the real propagation. 

Examples of the shape and evolution of stable and unstable solitons supported by the quintic-only (eo,i = 0) NL 
are displayed in Figs. [JJandlU In Fig. [JJ the comparison with the numerical findings demonstrates that the Gaussian 
ansatz fails to correctly predict the actual shape of the solitons in the case of the quintic nonlinearity, therefore only 
numerical results are reported in this section. For very narrow solitons, whose width is much smaller than the NL 
period [such as the soliton shown in Fig. UJb)], Eq. ^ is virtually tantamount to the equation with the constant 
nonlinearity coefficient, go = go + g\ = go + 1 (recall we fix g\ = 1), and, accordingly, the narrow solitons may be 
considered as Townes modes ([J), with go replaced by go- Therefore, the narrow solitons are unstable, see the central 
and right panels in Fig. [5J The instability manifests itself as a quick decay of the soliton (see the central panel in 
Fig- H]), or as the evolution towards the collapse, see the right panel in Fig. [5J Note the great difference in the 
scale of z in the panels of Fig. [2] demonstrating the stable and unstable evolution (i.e., the stable behavior persists 
indefinitely long, while the instability sets in quickly) . It is relevant to compare these scales with the diffraction length 
corresponding to the solitons displayed in Fig. [TJ which is Zdiffr ~ 0.1. On the other hand, it is relevant to stress that 
the GPE in the form of Eq. UJ may predict only the onset of the collapse, while its further development requires a 
more accurate (microscopic) description, in terms of the Hartree-Fock-Bogoliubov equations [401 ] . 

The numerically found dependence k(N) for the solitons in the quintic model is displayed in Fig. [3] for two different 
values of go. While these curves obey the VK criterion, direct simulations of the perturbed evolution of the solitons 
demonstrate that the solitons are stable only in narrow intervals of the values of N (recall that the VK criterion is 
only necessary, but not sufficient, for the stability). 

The numerically found stability region of solitons in the quintic model is located between the solid borders in Fig. 
21 for both positive and negative (self- focusing and defocusing) values of the constant quintic coefficient, go ^ 0. The 
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FIG. 3: The propagation constant vs. the norm for the soliton families in the quintic model (eo.i = 0) for go = (a) and 
go = 0.5 (b). The designated stability intervals were identified via direct simulations. 
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FIG. 4: Stability regions (located between the borders) for the solitons in the parameter plane of the quintic and CQ models 
(solid and dashed borders, respectively), shown separately for the self-defocusing (a) and self- focusing (b) sign of the constant 
part of the quintic coefficient, go- 



overall shape of the narrow region is quite close to 

N = (tt/2) V3/[2(1 + 5o )], (10) 

which is explained by the fact that the solitons are always found near values of N corresponding to the critical one, 
given by Eq. ([5]), with go substituted by the above-mentioned local value, go = go + 1. This argument also explains 
nonexistence of solitons at go < —1. The shrinkage of the stability stripe at large positive go is a manifestation of 
the relative weakness of the NL in this limit, and is consistent with the exponential weakness of the NL-induced 
stabilization predicted by the above analysis, see Eq. ([9]). Above the narrow stripe, the unstable solitons suffer the 
collapse, while below the stripe they decay, which may be realized as the delocalization transition, similar to that 
studied in systems with linear lattice potentials [HI, |4l| . 



IV. LOCALIZED MODES IN THE CQ (CUBIC-QUINTIC) MODEL 

Unlike the model with the quintic-only nonlinearity, the one with the CQ terms gives rise to broad solitons, as shown 
in Fig. [5] [see also Fig. (7{a) below]. Due to this fact, the predictions produced by the VA for the CQ model are in 
reasonable agreement with the numerical findings, as can be seen in Figs. [5] and [51 On the other hand, the agreement 
deteriorates at large values of N, when the quintic term dominates, driving the system to the regime considered in 
the previous section. 

Dependences k(N) in Fig. [6] satisfy the VK criterion, and simulations demonstrate that the solitons in the CQ 
model, in contrast to the quintic-only one, are always stable when they obey the criterion. Accordingly, the stability 
area found in the parameter plane of the CQ system, which is located between the dashed borders in Fig. 21 is much 
larger than its counterpart in the quintic model. The upper boundary of the area, which is determined by the collapse 
driven by the local quintic nonlinearity, is approximated by Eq. (fTTJj) . 

Another essential difference from the quintic-only model is that the stability area in the CQ system extends to 
go < — 1, where the quintic nonlinearity is completely defocusing, but the solitons can be supported by the cubic 
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FIG. 5: (Color online) Typical examples of broad stable solitons in the CQ model with go = 0.25, eo = ei = gi — 1: (a) 
k = 0.26, TV = 0.65; (b) k — 0.57, TV = 0.8; (c) k — 3.2, TV = 1.2. The dashed lines represent shapes predicted by the variational 
approximation for the same values of the parameters. 
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FIG. 6: Comparison of the variationally predicted and numerically found dependences k(N) in the CQ model, with eo = £i = 
g\ — 1, for go — (a) and go = 0.5 (b) . 



term. In this case, stable solitons develop side peaks, as shown in Fig. [71 and form odd-numbered stable in-phase 
bound states, see examples of three-, five- and seven-soliton complexes in Fig. [5] The transition from the higher-order 
soliton, featuring the side peaks, to the bound complexes occurs when the height of the additional peaks becomes 
equal to that of the central maximum. Bound states of three in-phase solitons were previously found in the ID 
cubic-NL model, but they were completely unstable [27 1. 

The present system also supports stable dipoles, i.e., bound states of two solitons with opposite signs, as well as 
stable dipolar complexes built of two multi-peak soliton clusterj^as shown in Fig. [5J The simplest two-peak dipoles 
are similar to their counterparts found in the cubic-NL model [27] . 




FIG. 7: Examples of a stable broad soliton (a) and stable solitons with symmetric side peaks (b, c) in the CQ model with 
e = e\ = gi = 1 and go = -1.05: (a) k = 2.88, TV = 19.5; (b) k = 5.22, TV = 21.3; (c) k = 9.23, TV = 22. 
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FIG. 8: Examples of stable in-phase bound states of three (a), five (b), and seven (c) solitons in the CQ model with eo = £1 = 
gi = 1, and g = -1.25: (a) k = 0.73, N = 20; (b) k = 1.12, TV" = 30; (c)k = 0.73, N = 40. 
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FIG. 9: Examples of stable dipole solitons (a, b), and of a stable cluster dipole (c), in the CQ model with eo = £i = 1- 
Parameters are: (a) g = 0.5, k = 1.1, N = 3.4; (b) g = -1.05, k = 2, N = 38; (c) g = -1.25, k = 1.1, N = 36. Note that 
the quintic term is self-defocusing in cases (b) and (c), but self- focusing, on the average, in (a). 



V. CONCLUSIONS 



This work aims to introduce the smooth sinusoidal NL (nonlinear lattice) acting on the critical quintic (x^) 
nonlinearity in the ID setting, as well as the system with the CQ (cubic-quintic) NL. The models approximate 
settings which are relevant to optics and BEC. Stability regions for solitons in these systems have been identified, 
being narrow in the quintic-only model, and much broader in the CQ one. In the latter case, the results may be 
predicted analytically by means of the VA (variational approximation), in the combination with the VK stability 
criterion. The narrow stability stripe in the quintic model, and its absence in the recently studied 2D models with 
sinusoidal NLs acting on the cubic nonlinearity (which is critical in the 2D geometry) were also qualitatively explained 
by means of an analytical approximation. A notable feature of the CQ NL is that it generates stable broad solitons, 
which may develop side peaks and form stable bound states, both in-phase and out-of-phase ones. 

It may be interesting to study other dynamical features in models of this type, such as mobility of the solitons. A 
relevant extension should also be consideration of two-component systems similar to the scalar models introduced in 
this work, cf. the analysis of two-component solitons in the cubic-NL system, which was reported in Ref. [42[ . 
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